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Abstract. This paper^ addresses the problem of scattering hne polarization and the 
Hanle effect in one-dimensional (ID), two-dimensional (2D) and three-dimensional (3D) 
media for the case of a two-level model atom without lower-level polarization and assum- 
ing complete frequency redistribution. The theoretical framework chosen for its formu- 
lation is the QED theory of Landi Degl'Innocenti (1983), which specifies the excitation 
state of the atoms in terms of the irreducible tensor components of the atomic density 
matrix. The self-consistent values of these density-matrix elements is to be determined by 
solving jointly the kinetic and radiative transfer equations for the Stokes parameters. We 
show how to achieve this by generalizing to Non-LTE polarization transfer the Jacobi- 
based ALT method of Olson et al. (1986) and the iterative schemes based on Gauss-Seidel 
iteration of Trujillo Bueno and Fabiani Bendicho (1995). These methods essentially main- 
tain the simplicity of the A— iteration method, but their convergence rate is extremely 
high. Finally, some ID and 2D model calculations are presented that illustrate the effect 
of horizontal atmospheric inhomogeneities on magnetic and non-magnetic resonance line 
polarization signals. 



1. Introduction 

The scattering line polarization, and its modification due to a weak mag- 
netic field — such that the Zeeman splitting is negligible compared with the 
line width (the so called Hanle effect; Hanle, 1924) — , sensitively depends 
on the anisotropy of the radiation field and on the magnetic field vector ge- 
ometry (Landi Degl'Innocenti, 1985; Stenflo, 1994). The solar atmospheric 
plasma is spatially inhomogeneous, with vertical and horizontal variations, 
not only in temperature, macroscopic velocity and density, but also in the 
orientation and intensity of the magnetic field (see Sanchez Almeida, 1999 
for new insights in this respect). It is thus clear that, in order to fully ex- 
ploit the Hanle effect as a diagnostic tool for weak magnetic fields, we also 



^Published in 1999 in the book Solar Polarization, edited by K.N. Nagendra & J.O. 
Stenflo. Kluwer Academic Publishers, 1999. (Astrophysics and Space Science Library ; 
Vol. 243), p. 143-156 
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need to address the problem of resonance line polarization and the Hanle 
effect in 2D and 3D media where the radiation field's anisotropy is different 
from that corresponding to the currently-assumed ID atmospheric models. 

To this end, this contribution begins presenting a formulation of reso- 
nance line polarization and the Hanle effect that we consider as the most 
suitable one for practical RT applications. It is based on the density-matrix 
theory for the generation and transfer of polarized radiation (see Landi 
Degl'Innocenti, 1983; 1984; 1985). In this paper we consider the standard 
case of a two-level model atom neglecting atomic polarization in its lower 
level (i.e. it is assTimed that the lower-level Zeeman sTiblcvels are equally 
populated and that there are no coherences among them). The quantities 
whose self- consistent values are to be determined are the irreducible tensor 
components of the density matrix {pg ) , which depend only on the spatial 
coordinates. The statistical equilibrium (SE) and RT equations to be solved 
are valid independently of whether we assume ID, 2D or 3D geometries. 

A summary of previous work done in the subject of the numerical so- 
lution of Non-LTE polarization transfer problems can be found in Trujillo 
Bueno and Manso Sainz (1999). In this respect, we should mention the 
recent work of Nagendra et al. (1998; see also their contribution in these 
proceedings) where the Hanle effect in ID is considered using a different 
theoretical approach. For information concerning the numerical solution of 
more general polarization transfer problems formulated with the density- 
matrix theory see Trujillo Bueno (1999). 

The outline of this paper is as follows. Section 2 presents the basic 
equations for the case of a two level atom without lower-level atomic po- 
larization and assuming complete frequency redistribution. Section 3 shows 
how the very efficient iterative methods of solution investigated by Trujillo 
Bueno and Fabiani Bendicho (1995) (Jacobi, Gauss-Seidel and successive 
over-relaxation) can be suitably generalized to the problem of Non-LTE 
polarized radiative transfer in the Hanle effect regime. Finally, in Sect. 4 
we present the results of some illustrative 2D Hanle-effect calculations for 
triplet lines and discuss the ensuing horizontal transfer effects. 

2. Basic Equations 

In scattering line polarization and Hanle effect problems, the quantum in- 
terferences (or coherences) between the magnetic sublevels of each atomic 
level must be considered. In order to properly take into account these ef- 
fects, we work within the framework of the polarization transfer theory 
based on the irreducible tensor components of the atomic density matrix 
(see Bommier and Sahal Brechot, 1978; Landi Degl'Innocenti, 1983). Since 
we are restricting ourselves in this paper to the standard case of a two-level 
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model atom neglecting lower-level atomic polarization only the K = Q = 
irreducible tensor component of the density matrix suffices to completely 
describe the lowcr-lcvcl excitation. For instance, if the lower level of total 
population A^; has angular momentum J; = one has that Pq{1) = Ni/^/S. 

However, with respect to the upper level we need to specify its total pop- 
ulation {pq{u)), its alignment {pq{u)), and two complex quantities (pi(n) 
and /92('u)) which take into account the coherences between the sub levels of 



the upper level. For example, for an upper level with = 1 

= -L(iVi + ATo + iV_i) = -^Nu, (1) 

p2(^) = -L(Ar,_2iVo + iV_i), (2) 

pUu) = -l=[pj^{l,0) - pjJO,-!)], (3) 

Pliu) = pj^il,-1), (4) 



where Ni {i = 1,0,-1) are the populations of the magnetic sublcvels 
M = 0, ±1 of the upper level, and pj{M, M') = <aJ M\p\a J M' > (with 
\aJM> the eigenvectors of the atomic Hamiltonian) are the elements of 
the density matrix in the standard representation (sec, e.g., Messiah, 1969). 
We point out that the orientation components {pq) of the density matrix 
are zero because we are assuming a static medium for which the radiation 
field that illuminates its boundaries has no circular polarization (see Landi 
Degl'Innoccnti et al, 1990). 

For each irreducible upper level tensor component pq with Q > 0, there 
exists another spherical component with Q < related to it through the 
conjugation property: 

where the symbol "*" means complex conjugation. We can thus choose the 
following linear combinations as independent variables: 

PQ =l[PQ + i-^fP-Q] = Re [p^] , Q > 0, (6) 

PQ = \^PQ- (-1)'^/^''q] = \PQ\ ' Q > 0' (7) 

where "z" is the imaginary unit. Finally, normalizing all these upper level 
unknowns to the total atomic population of the lower level, we find the six 
real unknowns of this problem: pQ, pg, pf , pf , p| and p| of the upper level. 

In the QED polarization transfer theory of Landi Degl'Innoccnti (1983, 
1984), the radiation field is described by its spherical tensor components: 

J^ = jdx^,j^h^ (8) 
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4 = j dx^,j^^[{^^^-l)hu + ^{^^^-l)QM (9) 

Jl = j dx ct>, j ^^e^x^l - + - iU,^] (10) 

J| = y'dx0,^^^e2^>^[i(l-/.2)4^^-^(l + /.2)Q,o-^/xt/,o] (11) 

where Ixn , QxCi and UxH are the Stokes parameters relative to the direction 
ri specified by the angles 9 and x defined as in Fig. (2b), jji — COS ^ and (px is 
the line profile, with x the frequency measured from the line center in units 
of the Doppler width. The physical meaning of these expressions is quite 
simple. Note that Jq is the well-known frequency integrated mean intensity. 
The other three irreducible tensor components of the radiation field are fre- 
quency and angular integrals of the three Stokes parameters, weighted by 
some angle-dependent quantities, and by the line profile (l)x- The radiation 
field tensor Jq measures the degree of vertical-horizontal anisotropy: it is 
positive when the radiation is predominantly vertical, negative if horizontal, 
and it vanishes at the bottom of the atmosphere where the radiation field 
is unpolarized and isotropic. Finally, the two other spherical tensor compo- 
nents are complex, and they measure the breaking of the axial symmetry of 
the radiation field through the azimuthal exponentials appearing inside the 
angular integrals. Therefore, they are zero in axially-symmctric media like 
ID plane-parallel atmospheres with a vertical magnetic field, or without 
any field at all. The Jq components with Q < Q can be obtained through 
a conjugation relation similar to the one stated in Eq. (5) for the atomic 
statistical tensors. With equivalent definitions to those of Eqs. (6)-(7), we 
obtain the six real quantities Jq, Jq, J\, jf , J^ and Jl. 

The SE equations that govern the six unknowns of this problem are: 



[1+ (5^(1-6)] 
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where = {2hvlJ^) [{2Ji + l)/^2J-^] p^. 

The Mjj-quantities of the magnetic operator M are coefficients that 

depend on the strength and orientation of the local magnetic field (see 

Table 1 in Landi Degl'Innocenti et al., 1990; for their explicit values), w^j^j^ 
is a numerical factor depending on the total angular momentum of the 
levels involved in the transition (see Table I in Landi Degl'Innocenti, 1984; 
and note that it is unity for a J; = and = 1 line transition), e is the 
coUisional destruction probability due to inelastic collisions, and B,^^^ is the 
Planck function. We point out that 6^^^ is the coUisional depolarizing rate 
due to elastic collisions measured in units of the Einstein A^i coefficient, 
with 5^^"^ = in the first equation. 

These SE equations have a clear physical meaning. The magnetic op- 
erator M couples the K = 2 statistical tensors among them. This is a 
local term because its Mjj coefficients only depend on the local value of the 
magnetic field. The second term in the r.h.s. of Eq. (12) is the radiative 
coupling term. It couples the atomic system with the radiation field and 
therefore it is highly non-local. Finally, the third term is the unpolarized 
thermal source. 

Since, as mentioned above, the orientation components are zero only 
three Stokes parameters are relevant in this problem: I^n, Qxn and Uxfi- 
Due to the fact that the lower level is assumed to be unpolarized, the 
radiative transfer equation for each Stokes parameter is decoupled from 
the others: 

-j — 2x^1 = Z^n — Sz, (13) 

dTx 

with Zj-n the Stokes parameter IxO.^ Qxii or UxSii and <iTx=—{xi4'x + Xc)(^s 
(with s the geometrical distance along the ray path and xi,c the line- 
integrated and continuum opacities). Although our code is very general, 
and takes into account the effect of a background continuum, for notational 
simplicity we will not consider it explicitly in the following equations. Thus, 
the line contributions to the source functions components Sz are: 

Sf^^ = + ^«£)_^^{-1=(3m2 - l)Si - V3n^/T^'{cosxSf - sinxSf) 



■2^/2 



+^{1 -fi^)icos2x Si -sin2x Si)}, (14) 



-^{1 + ^^^){cos2x Si -sm2x Si)}, (15) 
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aline 



W 



V3 { yr^(sin xSl + cos xSi) 
+//(sin2x5| + cos 2x^1)}, 



(16) 



where the S'q -quantities are given in terms of the density-matrix elements 
Pq as indicated by the expression given after Eq. (12). 

Since the transfer equations (13) are decoupled, it is straightforward 
to write the formal solution for IxUi QxCi and UxCi, as in the standard 
unpolarizcd case, through the monochromatic A^q operator (see Mihalas, 
1978). Thus, with the values of the density matrix elements we calculate Sj, 
Sq and Su, and then any suitable formal solution method of the standard 
RT equation can be used to calculate the Stokes parameters I^n, Qxii 
and Uxii at each grid point of the chosen spatial grid of NP points, for 
all the frequencies and directions of the chosen numerical quadrature. This 
allows us to write the radiation field tensors in the absence of a background 
continuum as: 



/ 



v 



TO 
Jf 

-P 

tin 



( Aoo 


Aoi 


A02 


Ao3 


Ao4 


Ao5 ^ 










/ 


-■-0 


\ 


Aio 


All 


A12 


Ai3 


Al4 


Ai5 












rr2 

-■-0 




A20 


A21 


A22 


A23 


A24 


A25 




s? 




+ 




j.2 




A30 


A31 


A32 


A33 


A34 


A35 










-tf 




A40 


A41 


A42 


A43 


A44 


A45 




Si 








rfi2 

^2 




\ A50 


A51 


A52 


A53 


A54 


A55 j 




Ui 


) 




\ 


rp2 
~ -"-2 


/ 



(17) 

where the NPxNP operators K^p are frequency and angular weighted av- 
erages of the standard K^q operator, and Jq , Sq and Tq are vectors of 
length NP, with the Tg components given by expressions similar to Eqs. 
(8)-(ll) but using the transmitted Stokes parameters due to the incident 
radiation at the boundaries instead of Ixn-, QxSi and Uxsi- The analytic ex- 
pressions of these A-like operators can be found in Manso Sainz & Trujillo 
Bucno (in preparation). As shown below, the only operator of relevance for 
our iterative approach is Aqo, which is nothing but the A operator of the 
standard unpolarized case: 



Aoo(i,j) = ^ J <l>x'^x jdO. Kuii,])- 



(18) 



It can be demonstrated that in plane-parallel atmospheres all the non- 
diagonal operators are zero except Aqi = Aio (see also Landi Degl'Innocenti 
et al, 1990; Nagendra, Frisch k Faurobert-Scholl, 1998), and that in 2D 
media A03, A05, A13, A15, A23, A25, A34, A45 and their symmetric ones 
are zero (Manso Sainz Sz Trujillo Bueno, in preparation). 
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3. Iterative methods of solution 

The most simple iterative scheme to solve this coupled set of equations 
is the A-iteration method. Given an estimate Sq of the unknowns (i.e. 
of the Sq tensors at all the spatial grid-points) solve formally the transfer 
equations to calculate the corresponding six Jq"'''^ radiation field tensors at 
each spatial grid-point "z". Then, introduce these JQ°^'^{i) values into the 
SE equations at each spatial grid-point independently and get an improved 
new set of 5q -values. As shown by Trujillo Bueno and Manso Sainz (1999), 
this A-iteration method can be used as a reliable solution method if one 
initializes with the self-consistent pg-valucs corresponding to the case in 
which polarization phenomena are neglected. 

In order to develop iterative methods characterized by an extremely high 
convcu'gcncc rate, and that can be applied to find the self-consistent solution 
independently of the chosen initialization, we need to account implicitly for 
some of the "new" values of the unknowns pg-elements. To this end, in 
the following we generalize to the Hanle effect regime the Jacobi-based ALI 
method of Olson et al. (1986) and the iterative methods of Trujillo Bueno 
and Fabiani Bendicho (1995) that are based on Gauss-Seidel iteration (see 
also Trujillo Bueno and Manso Sainz, 1999). 

In order to derive a Jacobi-based ALI scheme we do the same as with the 
A-iteration method, except that in order to calculate Jg (?) we use the new 
value of Sq (instead of the old one) at the grid point "i" being considered. 
Since this "ne?i/' value is not yet known, we implicitly write it only in the 
expression of Jq, i.e. we write 

jO(z)«Jo°°''^(i) + Aoo(i,?)<550(i), 

4{:i)^J^"^'{i), 

where (55'g(i) = S'g"'^"'(i) — S'o"''^(i). After substitution of these expressions 
for the Jq (i) quantities into the SE equations, we obtain at each grid-point 
"i" independently a system of six equations with six unknowns that can be 
solved easily to find the new values of the six statistical tensors Sq (i). We 
point out that this is equivalent to applying the operator splitting technique 
to the Ago operator only. It can be demonstrated that no gain is obtained 
if the splitting is apphed to the whole set of 36 operators of Eq. (17) (see 
next subsection below). 

We emphasize that, at each iterative step, the SSq corrections are made 
point by point. Thus, a better idea than Jacobi's method would be to apply 
the method based on Gauss-Seidel (GS) iteration of Trujillo Bueno and 
Fabiani Bendicho (1995), i.e. to do the same as with the A-iteration method, 
except that in order to calculate Jq (i) we use the new values of Sq (instead 
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of the old ones) at the grid point "i" being considered and also at the grid- 
points l,2,...,z — 1 that have been previously considered as we advance from 
the atmosphere's boundary at which the first correction was made. Since 
5*0 "'^ (i) is not yet known, we imphcitly write it only in the expression of 
Jq (i), i.e. we write 

J0(i) « jOold^ndnen.^^ ^ Aqo (i, (i) , 

4(^)«J2o^^(^), (20) 

where jOoidandneiiJ^-^-j -g ^^iq avcragcd mean intensity calculated using the 
^^new" values of Sq at the grid-points l,2,...,i — 1 and the "oW S'g -values 
at points i,i + l,i + 2, NP, with NP the total number of points of the 
spatial grid (see Trujillo Bueno and Fabiani Bendicho, 1995: for details). 

After substitution of these expressions for the Jq (i) quantities into the 
SE equations, we obtain at each grid-point "i" a system of six equations 
with six unknowns that can be easily solved to find the new values of the 
statistical tensors SQ{i). By implementing this GS-bascd iterative scheme 
as suggested in the conclusions of the paper by Trujillo Bueno and Fabiani 
Bendicho (1995) the total computational work required to achieve the self- 
consistent solution is a factor 4 smaller than with the previous Jacobi-based 
method. 

An extra important improvement of the convergence rate can be achieved 
by multiplying each 6Sq GS-correction by a numerical factor (a;) lying be- 
tween 1 and 2. This is the successive over-relaxation (SOR) method, i.e. 

Ssf""^ = udS^""' (21) 

The optimal value of u can be easily found (see Trujillo Bueno &: Fabiani 
Bendicho 1995); however, a good convergence rate can be obtained by just 
choosing uj = 1.5. 

Figure 1 shows the convergence rates of the six Sq unknowns, i.e. it 
shows the variation with the iteration number of the maximum relative 
change (Re) in these quantities. The problem considered was the Hanle 
effect in a 2D model atmosphere with e = 10~^, S^'^'^ = and with a non- 
vertical magnetic field vector. The solid lines refer to the Jacobi-based ALI 
method. We point out that the Sq tensor elements converge at the same 
rate as Sq. The reason for this is that all the Jq radiation field tensors 
are dominated by the Stokes I parameter, and this specific intensity is 
basically set by the value of Sq (see Trujillo Bueno and Manso Sainz, 1999). 
Furthermore, the convergence rate of Sq is the same in polarized (with 
or without magnetic fields) and in unpolarized problems as shown by the 
dotted line in Fig. 1 that cannot be distinguished from the Sq solid-line (see 
also Faurobert-Scholl et al., 1997). It is important to point out that this is 
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Figure 1. The convergence rate of the six Sq unknowns for the Jacobi method (soUd 
hnes), and of Sq for the Gauss-Seidel and the SOR method with u> = 1.5 (dashed hnes). 
See the text for the meaning of the dotted hne. 

a general behaviour and it is independent of the geometry of the medium 
(ID, 2D or 3D). 

The convergence rates for the Gauss-Seidel and SOR iterative methods 
are also plotted in Fig.l. For clarity reasons the convergence rates for the 
Sq elements have been omitted because they present a similar behaviour. 
As seen in the figure OTir GS method is four times faster than Jacobi, while 
our SOR method would be a factor 10 if the calculation had been performed 
with the optimal a;- value (see also Trujillo Bueno and Manso Sainz, 1999). 

4. Some illustrative examples 

In the following we show some illustrative examples of the Hanle effect in 
2D atmospheric models comparing the results with the corresponding ID 
case. Although our code is very general and can deal with realistic 2D and 
3D scenarios with horizontal periodic boundary conditions, here we will 
restrict ourselves to sinusoidal fluctuations of the Planck function along 
the horizontal X-axis: 



with k^. = 2'k/L the horizontal wavenumber, and L the horizontal wave- 
length of the thermal inhomogeneities. All the geometrical distances are 
measured in units of the opacity scale height (7^^ ~ 100 km) . We assume 
a gaussian line absorption profile and consider different A;a;-values, being 
kx = Q the plane-parallel ID limit. 

Figure 2 is presented to facilitate the understanding of the geometry of 
the problem and to indicate the chosen positive and negative directions of 



Bi, = Bu{z) + ABjy cos{kxx) 



(22) 
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Figure 2. (a) Schematic geometry of a 2D atmosphere with a sinusoidally varying 
temperature inhomogeneity. (b) Show the angles defining the Une-of-sight (ft) and the 
polarization unit vectors ei and e2. The positive QxSi direction is along ei. (c) Angles 
defining the magnetic field orientation. 

the Stokes Q-parameter. In Fig. 2a the shaded and white "slabs" simply 
aim at visualizing the coolest and hottest regions of the assumed 2D model 
atmosphere, respectively. Note, however, that the chosen 2D medium is not 
composed of such embedded slabs, since it is characterized by the horizon- 
tal temperature fluctuations given by the previous equation. These "slabs" 
should be understood as infinite along the Y-direction. Fig. 2b shows the 
angles that specify the direction of propagation of the ray under consider- 
ation, while Fig. 2c gives the angles that determine the orientation of the 
magnetic field vector. As seen in Fig. 2a the Y-Z plane at X=0 is dividing 
the hottest "slabs" in two equal halves. We will show the emergent polar- 
ization profiles for simulated observations made along the hottest "slabs" 
(i.e. at the horizontal position X=0 and for a line of sight with x = 90"). 
Thus, the positive direction of the Stokes Q-parameter lies along the slabs, 
and the negative one is perpendicular to them. 

In Fig. 3 we consider the case of resonance line polarization in this 2D 
atmosphere. When observing at disk center (i.e. at /x = cos^ = 1) we find 
that there is no polarization signal corresponding to the plane-parallel ID 
limit (i.e. for = 0), as it must indeed be the case because the radiation 
field of an unmagnetized ID medium has axial symmetry, and the radia- 
tion field tensors Jf and j| are then zero. As we increase the horizontal 
wavenumber kx (i.e. the parameter that measures the degree of horizontal 
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Figuixi 3. Emergent Q/I profiles at disk center (/u = 1) for different liorizontal wavenum- 
ber fcx-values. Here AB^/B^ = 0.1. 

inhomogeneity of our 2D model) the radiation field looses its axial symme- 
try, and the polarization signal Q/I starts to increase accordingly. However, 
when the horizontal temperature inhomogcneitics are smaller than about 
2000 km (i.e. for wavenumbers kx > 0.3), the polarization signal decreases 
as the wavenumber kx further increases. This is because we start then to 
approach the limiting case of an atmosphere composed of optically thin 
irregularities, for which the radiation field recovers the axial symmetry 
characteristic of a ID medium. 

Fig. 4 shows the emergent Q/I and U/I profiles for simulated high- 
spatial resolution observations with a line of sight having x = 90" and made 
close to the solar limb (/U = 0.1) in a 2D atmosphere with /S.Bi,/By = 0.1 
and kx = 0.3 (i.e. with Lfa 2000 km). The magnetic field direction is defined 
through the azimuthal and polar angles and 9b (see Fig. 2c), and its 
strength through the parameter F = 0.88gjB/Aui, with gj the upper level 
Lande factor, B the magnetic field in Gauss, and A^i the spontaneous 
emission Einstein coefficient measured in units of lO^s^^. Since the line- 
of-sight lies along the slabs (i.e. x = 90°)) we find that in the absence 
of magnetic field (i.e. for F = 0) U/I is zero due to symmetry reasons. 
The two uppermost panels show how a magnetic field parallel to the slab 
(i.e. with (j)B = 9b = 90°) substantially changes the emergent polarization 
signal. However, when the magnetic field is perpendicular to the slab (see 
the lowermost panels with 9b = 90° and (j)B = 0^), the change with F 
of Q/I is the smallest one and U/I always remains zero. The two central 
panels show the intermediate case with (f)B = 45°. The dashed lines show 
the results for the plane-parallel kx = case. 

Fig. 5 shows the emergent fractional linear polarization Q// at /x = 0.1 
and X = 90°, when the Planck function varies sinusoidally (with kx = 0.3 
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Figure 4- Emergent Q/I and U/I profiles in the assumed 2D atmosphere (solid lines) 
with wavenumber kx = 0.3 and AB^/Bi, = 0.1. The line of sight observation is at /n = 0.1 
and X = 90°. Dashed lines indicate the ID case. 



and AB^/B^ = 0.2), and there exists additionally a horizontal fluctuation 
in the line opacity xi that is anticorrelated with that of the Planck function, 
i.e. 

Xi = Xi[l + "cos(fca:x)], (23) 

with kx = 0.3 and a = —0.2. The curves situated on the left-hand-side of 
the figure show the ID and 2D results corresponding to the zero magnetic 
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Figure 5. As Fig. 4 but taking AB^/Bv = 0.2 and including the horizontal opacity 
inhomogeneities given by Eq. (23). The magnetic field has F = 1,6b = 90°, and 4) = Q° . 
Solid lines show the 2D case, and dashed-lines the ID one. 



field case, while the curves on the r.h.s. correspond to a case of a magnetic 
field with T = 1, Ob = 90°, and 4)b = 0°. For this magnetic and line- 
of-sight geometry U/I = 0. It can be seen that, for this particular model 
atmosphere and geometry, a Hanle-effect diagnostic of very high spatial 
resolution observations would lead to an underestimation the magnetic field 
strength for interpretations based on the plane-parallel ID approximation. 

5. Conclusions 

We have developed a Hanle effect code that allows the numerical simulation 
of resonance line polarization signals in the presence of weak magnetic fields 
in ID, 2D and 3D media. The governing equations have been formulated 
working within the framework of the density matrix polarization transfer 
theory of Landi Degl'Innocenti (1983, 1985). These SE and RT equations 
are the same independently of whether we are considering ID, 2D or 3D 
atmospheric models. The six pg-unknowns of the problem are neither fre- 
quency nor angle dependent, since they only vary with the spatial position. 

Three different iterative schemes that were originally developed for RT 
applications in the unpolarized case have been generalized to solve this 
set of equations: Jacobi (ALI), Gauss-Seidel and SOR (Trujillo Bueno & 
Fabiani Bendicho, 1995, Trujillo Bueno & Manso Sainz 1999). This kind of 
iterative methods does not make use of any matrix inversion, and essentially 
maintain the A- iteration simplicity. The only difference between the ID, 
2D and 3D versions of our Hanle effect code lies in the formal solution 
routine that calculates the radiation field tensors from the current values 
of the density-matrix elements. To this end, in 2D we use the formal solver 
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developed by Auer, Fabiani Bendicho and Trujillo Bueno (1994) and in 3D 
we use the one presented at this workshop by Fabiani Bendicho and Trujillo 
Bueno (1999). 

We have also shown some Hanle effect results for ID and 2D media. 
These calculations illustrate how weak magnetic fields and horizontal ra- 
diative transfer effects compete to modify the scattering line polarization 
signals expected from plane-parallel ID atmospheres. Thus, further careful 
investigations must be done in order to separate both effects, when diag- 
nosing weak solar magnetic fields via the Hanle effect. This type of future 
studies should be done thinking in the interpretation of low spatial reso- 
lution scattering line polarization observations. Our numerical approach is 
very efficient and suitable to investigate scattering polarization signals for a 
variety of atmospheric models having any desired temperature, density and 
magnetic field vector variations. Another useful research that can be done 
with our Hanle effect codes concerns the simulation of polarization signals 
emerging from realistic MHD and semi-empirical 2D and 3D models. 
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